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ABSTRACT 

We have calculated Fe II emission line strengths for Active Galactic Nuclei 
Broad-Line Regions using precise radiative transfer and Iron Project atomic 
data. We improve the treatment of all previously considered excitation 
mechanisms for the Fe II emission, continuum fluorescence, collisional excitation, 
fluorescence by self-overlap among the iron lines, and fluorescent excitation by 
Lya. We demonstrate that Ly« fluorescence is of fundamental importance in 
determining the strength of the Fe II emission. In addition to enhancing the 
ultraviolet and optical Fe II flux, Lya fluorescence also results in significant 
near-infrared Fe II emission in the AA8500 - 9500 A wavelength range. New 
observations are suggested to probe this effect in strong Fe II emitting quasars. 

Subject headings: quasars: emission lines — radiative transfer — line: formation 
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1. Introduction 

One of the most puzzling instances of Fe II line formation occurs in the spectra 
of Active Galactic Nuclei (AGN) with Broad-Line Regions (BLR) which exhibit strong 
ultraviolet Fe II emission (UV = AA2000 - 3000 A) and often strong optical Fe II emission 
(OP = AA3000 - 6000 A) as well. The widths of the iron lines indicate formation in the 
dense BLR, yet standard photoionization calculations, successful in interpreting many 
other BLR features, fail to account for the strength of the Fe II emission (Wills, Netzer & 
Wills 1985; Joly 1993). Observed total Fe II/H/3 ratios are 12, whereas standard models 
predict ratios less than 8, typically 3-5, with the largest values resulting from enhanced 
iron abundances. There is also a class of "superstrong" Fe II emitters with Fe II/H/3 « 30 
which seem unable to be explained by traditional photoionization modeling (Joly 1993; 
Graham et al 1996). 

Excitation mechanisms advanced for the Fe II emission include continuum fluorescence 
via the UV resonance lines, self-fluorescence via overlapping Fe II transitions, and collisional 
excitation, with the latter thought to contribute the bulk of the Fe II emission. However, 
Wills et al (1985) demonstrated that model BLR calculations including these excitation 
mechanisms cannot account for the observed strengths of the Fe II lines. Elitzur & Netzer 
(1985, hereafter EN) extended this calculation to include Lya fluorescent excitation but 
found its contribution to be negligible. Nevertheless, Penston (1987) argued that strong 
indirect observational evidence for the importance of Lja fluorescence is found in the 
presence of unexpected UV Fe II multiplets in the spectrum of the symbiotic star RR Tel 
which seem attributable only to cascades from higher levels pumped by Lya fluorescence. 
The emission nebulae of symbiotics offer densities and ionization parameters similar to 
the BLRs of AGN. Recently, Graham et al (1996) have identified emission from these UV 
Fe II multiplets expected to be preferentially strengthened by Lja in the spectrum of the 
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ultra-strong Fe II emitter 2226-3905. 

In this Letter, we improve the modeling of all of the mentioned excitation mechanisms 
by using new, accurate atomic data and precise radiative transfer, and we demonstrate 
that fluorescent excitation by Lya can be of prime importance. A key feature of Lya 
fluorescent excitation of Fe II is significant near-IR Fe II emission in the wavelength range 
AA8500 — 9500 A. We suggest that a currently well-known, but unidentified, emission feature 
near A9200 in AGN spectra may be, at least in part, a spectral signature of this process. 

2. Calculations 

Fe I - IV were considered with Fe I, Fe III, and Fe IV each represented by a ground state 
photoionization cross section and total, unified radiative plus dielectronic recombination 
rate from the Iron Project calculations (see Bautista & Pradhan 1998). The Fe II atomic 
model consists of 262 fine-structure levels with energies from Johansson (1978; 1992 private 
communication). The model atom is complete in quartet and sextet levels for energies 
below 0.849 Rydbergs, allowing full inclusion of all Lya fluorescently excited transitions 
from a 4 D. This atom is a subset of the much larger 827 level Fe II atom used by Sigut & 
Pradhan (1996). 

Transition probabilities for Fe II were drawn from 3 sources, Fuhr, Martin & Wiese 
(1988), Nahar (1995), and Kurucz (1991, private communication). Over 3,400 Fe II 
radiative transitions were included in the non-LTE solution. Photoionization cross sections 
for all Fe II levels were taken from Nahar & Pradhan (1994). A level-specific radiative 
plus dielectronic recombination coefficient was used for each Fe II quartet and sextet level. 
For collisional excitation of Fe II, we adopted the R-matrix results of Zhang & Pradhan 
(1995) and Bautista & Pradhan (1998) and used the Gaunt-factor formula for all remaining 
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high-lying, dipole transitions. Charge transfer recombination rates to Fe II and ill were 
taken from Kingdon & Ferland (1996). Charge transfer ionization of Fe II was also included 
as Fe III + recombination is to the Fe II ground state (Neufeld & Dalgarno 1987). 

The coupled equations of radiative transfer and statistical equilibrium were solved 
with the accelerated lambda-iteration method of Rybicki & Hummer (1992). As overlap 
between the Fe II transitions is complex, the full preconditioning strategy was implemented. 
Complete frequency redistribution (CRD) over depth- dependent Doppler profiles was 
assumed for all iron transitions. The calculations were performed in plane-parallel, one 
dimensional, continuum illuminated clouds, each computed for a given ionization parameter 
and total particle density, assuming constant total pressure. The run of T e and N e , the 
photoionizing radiation field, and the continuous opacities and source functions were 
obtained with the cloudy code (Ferland 1991). This cloud structure was fixed during the 
iron non-LTE calculation. The solar iron abundance was used, \og(N Fe /Nii) + 12 = 7.477. 

The Lja source function can be written without significant error as 



because stimulated emission is not important. The quantity in brackets is the usual 
CRD source function. The absorption profile, <f) v , was taken to be a depth-dependent 
Voigt profile with a width set by natural broadening. The emission profile, tp u , however, 
depends on the radiation field (Mihalas 1978). We have retained the cloudy H I Is 
and 2p level populations, but have explicitly computed the emission profile in the partial 
coherent scattering approximation as formulated by Kneer (1975). Assuming a fraction 7 C 
of Lya photons are coherently scattered in the atom's rest frame, the ratio of emission to 
absorption profile is 




(1) 
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where J v is the mean intensity, and symmetric function that effects the transition 

from CRD in the line core to coherent scattering in the wings. We have taken the form 
of a v i }V from Vernazza, Avrett & Loeser (1981) and adopted 7 C = 0.98. We solved the 
radiative transfer problem implied by equation @) separately with an iterative technique, 
and then the Lya source function and opacity were incorporated into the appropriate 
monochromatic source functions during the solution of the radiative transfer equations 
for the Fe II transitions. We have verified our solution of equation (^) by matching the 
Lya profiles tabulated for BLR clouds by Avrett & Loeser (1988) using their models and 
hydrogen populations. 

3. Results 

Figure [l] compares the predicted Fe II flux with and without Lya pumping in BLR 
models with two different ionization parameters (Ui on ), 10~ 2 and 10~ 3 . These values are 
chosen to be representative of typical Fe II-emitting clouds, but do not exhaust the range 
of possibilities. Both models assumed a total particle density 4 • 10 9 cm -3 , a total column 
density of 10 23 cm -2 , and zero microturbulence. The shape of the photoionizing continuum 
was that of Mathews & Ferland (1987). 

For [/i on = 10~ 2 , the influence of Lya is large, more than doubling the total emitted 
flux, giving Fe II/LL3 = 2.4 for UV emission and 0.78 for optical emission, as compared to 
1.1 and 0.29, respectively, without Lyct (the H/3 flux was taken from the cloudy model). 
The strongest flux enhancements occur as noted by Penston (1987). The most important 
fluorescent excitation rates are from low-lying a 4 D e levels to the 4 D°, 4 F°, 4 P°, 4 S°, and 
6 F° symmetries of the 3d 6 ( 5 D)5p configuration (Johansson & Jordan 1984). Cascades from 
these levels to e 4 D and e 6 D produce strong emission lines in the range AA8500 — 9500 A 
(see Fig. |2|). In turn, cascades from e 4 D and e 6 D to odd parity levels near 5 eV enhance 
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multiplets 399, 391, 380 (~ A2850), 373 (~ A2770), and 363 (~ A2530); these lines are the 
unexpected UV multiplets in the RR Tel spectrum noted by Penston and the multiplets 
identified by Graham et al (1996) in 2226-3905. Penston estimated that « 20% of the UV 
Fe II flux comes out in these multiplets, which compares favorably with the ~ 30% found 
in the current model. Subsequent cascades then enhance the optical and UV fluxes from 
transitions among the lower levels. While transitions from a 4 D within ±3 A of Lya are the 
most important in the fluorescent process (see the listing of Johansson & Jordan 1984), 
our treatment of Lya includes the entire profile within ~ ±50 A. Including pumping only 
within ±3 A reduces the flux ratios to 2.1 for UV emission and 0.70 for optical emission for 
the case of U ion = 10~ 2 . 

For [/ion = 10~ 3 , the effect of Lya pumping is not as large, particularly on the UV 
flux. The Fe II/H/3 ratio generally increases with decreasing ionization parameter (even 
without Lya pumping) because of the larger sensitivity of the hydrogen level populations 
to temperature. 

As noted by Penston (1987), and confirmed by the current calculation (Fig. 0), strong 
observational evidence for Lya pumping of Fe II would come from the identification of the 
initial decays from the pumped 5p levels to e 4 D and e 6 D4 which give rise to transitions in 
the AA8500 — 9500 A wavelength range. We note that there is a prominent unidentified 
feature near A9200 in quasar spectra as discussed by Morris & Ward (1989) and Osterbrock, 
Shaw & Veilleux (1990). An identification as solely H I Pa9 A9229 would render it 
anomalously strong compared to Pa8 and PalO. Morris & Ward tentatively assign this 
feature to Mg II 4p 2 P° /2 1/2 - 4s 2 Si /2 (A9218, 9244) as formed by cascades from 5p 2 P^ /2A/2 
pumped by Ly/3. They discuss possible contributions from Fe II but discount this possibility 
based on the apparent absence of a line at A8927. Our calculations indicate the A8927 
Fe II line is significantly weaker than the main feature near A9200 (see Fig. Q). We have 
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performed a preliminary non-LTE calculation for Mg II including Ly/3 fluorescence in a 
manner similar to that described for Lya (but with 7 C = 0.4) and find that while the Mg II 
transitions in this region are strengthened by Ly/3 by nearly a factor of 5, they are still 
only minor contributors to the flux in the A9200 region (< 5%). Hence, the A9200 feature 
may be, in part, a spectral signature of Lya pumping of Fe II. As noted by Osterbrock et 
al, observations in this wavelength region are hampered by telluric H 2 absorption bands, 
which require careful removal, and, for available observations, falling CCD responses in the 
near-IR. New observations of this wavelength region are warranted. 

The current efficiency of Lya pumping is in disagreement with the EN calculation. 
An important difference is in the adopted oscillator strengths of the pumped transitions 
from a 4 D. The Kurucz (1981) values are roughly a factor of 10 smaller than those obtained 
by either Nahar (1995) or Kurucz (1991, private communication). Reducing the oscillator 
strengths for the pumped transitions within ±3 A of Lya by a factor of 10 results in a 
strong reduction in the effect of Lya fluorescence, giving for the U{ on = 10~ 2 model of 
Figure [l|, Fe II/H/? ratios of 1.28 for UV emission and 0.36 for optical emission. The near-IR 
emission (Fig. |2|) is also substantially reduced from Fe II/H/5 = 0.326 to 0.061. EN also 
used a simple escape probability treatment of line overlap, and a simpler treatment of Lya, 
so significant differences with the current work are perhaps expected. Finally, we note that 
our atom is essentially complete in the large number of quartet and sextet Fe II energy 
levels below the 3d 6 ( 5 D)5p configuration, and that our branching ratios from the pumped 
levels are expected to be considerably more accurate than those used by EN. 

4. Conclusions 

Using precise radiative transfer and accurate atomic data for Fe II, we demonstrate 
that Lya fluorescent excitation can more than double the Fe II flux in the UV and optical 
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regions. We also predict Fe II emission in the AA8500 - 9500A region as a direct result of 
Lya fluorescence. Both of these results can be viewed as theoretical confirmation of the 
role of Lya fluorescence suggested by Penston (1987) on the basis of observations. New 
observations of the AA8500 — 8500 A wavelength range in strong Fe II emitters are highly 
desirable, and would provide valuable constraints on modeling Lya fluorescence. Potentially 
new information on the Lya profile within single BLR clouds could be obtained, something 
which is impossible based on observations of the highly broadened Lya profile itself. 

A detailed discussion of the variation of the Fe II emission with the physical parameters 
defining the BLR clouds (using a much larger 827 level Fe II atom which includes the 
doublet spin system), and the implications for the 'Fe Il/H/3' problem, will be presented by 
Sigut & Pradhan (1998). 

We are also extending our calculations to overcome some of the limitations of the 
present work, in particular, to incorporate all important ions into a full solution of the 
non-LTE radiative transfer problem for the thermal structure of AGN BLR clouds. This 
will enable the Lya-Fe II problem to be solved within a self-consistent framework. 

This work was supported by the US NSF through grants to the NCAR High Altitude 
Observatory and to AKP (PHY-9421898). This work was begun while TAAS was an 
NSERC postdoctoral fellow at HAO. TAAS wishes to thank J. M. Marlborough and J. D. 
Landstreet for support through their NSERC grants. 
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Fig. I. — The predicted Fe II spectrum. The spectra were convolved with a Gaussian of 1/e 
width of 2000 kms" 1 . Covering fractions of the illuminated and shielded cloud faces were 
assumed to be 5% each. The ratios of Fe II fluxes with H/5 are computed with the optical 
(OP) and ultraviolet (UV) wavelength regions defined in the introduction. 

Fig. 2. — An enlargement of the near-IR Fe II spectrum for the U\ on = 10 -2 model of 
Figure 1. The spectra were convolved with a Gaussian of 1/e width of 500 kms -1 . 'Lya 
(1/10)' indicates that the Aji values of the pumped Fe II transitions within ±3 A of Lya 
were scaled by 0.1. Fe II(IR) refers to the flux between 0.7 — 1.0 /zm. 
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